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Abstract. This article reports on the results of simulations conducted to assess the performance of a modal Multi-Conjugate 
Adaptive Optics (MCAO) system on a 10m telescope with one Deformable Mirror (DM) conjugated to the telescope pupil and 
a second DM conjugated at a certain altitude above the pupil. The main goal of these simulations is to study the dependence of 
MCAO performance upon the altitude of the high-altitude conjugated DM and thereby determine its optimal conjugation. The 
performance is also studied with respect to the geometry of the Guide Star constellation when using constellations of Natural 
Guide Stars (NGS), which are rare, or constellations of Laser Guide Stars (LGS) which would allow large sky coverage. 
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1. Introduction 



Adaptive Optics (AO) was first proposed by Babcock ("19531) 
as a means to compensate for the loss of spatial resolution of 
telescopes due to the distu rbing effects caused by atmospheric 
turbulence ('Roddiei"l98l') and thereby obtain diffraction lim- 
ited resolution. A major limitation of AO is the small size of 
the Field of View (FoV) over which the correction is effec- 
tive, referred to as the isoplanatic patch. Under good seeing 
conditions a typical value for the isoplanatic radius is just a 
few (2-3) arcsec in the visible (Q.Sfim) and about 12-18 arc- 
sec in the near infrared (2.2fj.m). The concept of MCAO as 
a technique t o deliver large corrected FoVs w as pro posed by 
IPickd ill 9751) and later developed by BeckersI lll988l) with the 
main goal of increasing the corrected FoV. It did not however 
receive much attention at that time as the AO community fo- 
cused on conceptual problems such as obtaining adequate sky 
coverage with and without the use of Laser Guide Stars (LGS) 
(Sandler et al. 1994; Le Louarn et al. 1993) ™d the problems 
associated with LGS-based AO systems, in particular tip-tilt 
indetermination and the cone effect. During the early '90s the 
concept of atmospheric turbulence reconstruction in 3-D (i.e. 
tomography) was proposed as a solution to overcome the prob- 
lems associated with sky coverage and the cone effect in LGS- 
based AO systems (Tallon & Foy 1990). Shortly afterwards, 
the theoretical foundations of MCAO wer e developed (see e.g. 
IJohnston & W elsh'l994VEllerbr oekll994i) . In the last few years 
it has been recognised that tomography and MCAO offer the 
possibility of attaining large corrected FoVs and there is a large 



on-going effort in MCAO 


on both the theoretical side (see 




^^agazzoni et al. 1999; Fusco et al. 2000; Ragazzoni et alj 


Eoof 


i iRigaut et alJ 


l200(i iFuscoetalJ 2001) and using nu- 


merical simulation (see e.g. 'Flicker et al. 2000; Carbillet et alj 



2001a; Conan etal. 2001; Bello et al. 2001; Ellerbroek 2(m: 
Tordi et alJl200ll: IXo kovinin et al. 200 It iFemenia et alJl2()0lt 
Ellerbroek & Rigaut .200kiLe Louarni<2002h to assess the per- 
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formance of MCAO systems. 

Conventional AO systems employ a single Deformable 
Mirror (DM) usually conjugated to the telescope entrance 
pupil. MCAO requires at least a second DM conjugated to some 
altitude above the pupil. The addition of a second DM repre- 
sents a modest increase in complexity (of course the wave-front 
sensor also has to be modified). Such a dual-conjugate mode is 
considered as a short term upgrade to the AO system bei ng de- 
signed for the 10m GTC telescope ( iDevanev et alJ20o3) at the 
Observatorio del Roque de los Muchachos (ORM, Spain). In 
this work we present numerical simulations of a dual-miiTor 
modal MCAO system for the case of a 10-m telescope. The 
main goal of our simulations is to determine the optimal conju- 
gation altitude of the high-altitude DM (from now on referred 
to as DM2) assuming that the low-altitude DM {DM\) is con- 
jugated to ground level. Our numerical simulations investigate 
the dependence of the optimal high-altitude conjugation alti- 
tude on the nature of the GS constellation (either NGS or NGS 
plus LGS constellations are considered), the geometry and the 
number of reference sources in the GS constellation. These 
simulation results will be com |)ared with theoretical pre dic- 
tions based on the estimates of 'Tokov inin & ViardI 1I2OOII) for 
the optimal conjugation altitude of DM2. 
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2. MCAO model 

2.1. Modal tomography approach 

In modal tomography (iLawrence &Chowl 11984 
[Raaazzoni et al. 1999) one seeks to estimate the coefficients of 
Zernike expansions of the wavefronts over metapupils defined 
at the conjugate altitudes of the DMs. The metapupils are 
circles enclosing the footprints of the beams from the guide 
stars. If /7(g) denotes the Zernike polynomial decomposition 
of the incoming wavefront when looking at guide star g in 
the direction "^(g) (i.e. W(^; ~§ (^)) = 2; biXg)Zi(^ jR), with R 
the pupil radius) then the Zernike polynomial decomposition 
of the wavefront defined on each of the metapupils to be 
reconstructed at heights hi (I - I . . .L) is |"^'''}^_^ ^ which is 

related to [^(g)]^^^ ^ by : 

"^^ = p-s'\ (1) 

where the subscript ^ denotes transpose and the difi'erent terms 
in Eq. Q are: 



b of the GS wavefronts are retrieved from real data. In a real 
wavefront sensor the data will be wavefront local slope mea- 
surements (as in the case of Shack-Hartmann and pyramid 
wavefront sensors) or wavefront local curvature values (as in 
the case of the curvature sensor). We require a reconstruction 
process from real measurements to the metapupil Zernike co- 
efficients 7? without having to first reconstruct the GS Zernike 
coefficients b. This is achieved by means of the interaction 
matrix M of the MCAO system linking wavefront sensor out- 
puts fi?) and metapupil Zernike coefficients: "i? - M"^. While 
the P matrix in Sect. 12.11 is obtained from purely geometric 
considerations, building the interaction matrix M requires ex- 
plicit knowledge of the WFS such as the extended size of the 
GS spot as seen by the WFS, pixelization effects in the WFS or 
the fact that the WFS sees only a limited FoV. This leads us to 
use interaction matrices generated by the simulation in which 
we sequentially place different Zernike polynomials on each 
layer and register the local tilts produced on each subaperture 
of each Shack-Hartmann Sensor (SHS) in our MCAO system. 
In the simulations using LGSs we consider that there is a global 
tilt subtraction for the wavefront coming from each LGS. 



= t(2)("^(2)) ... t(G)C^(G))] 

P('^(,);/2i) P("^(i);/Z2) ■■■ P(t(,);/2L) 
P("^(2);/!i) P('^(2);/!2) ■■■ P("^(2);/!l) 

P('^(G);/ii) V{d(G)\h2) ■■■ P0(Gy,hL) 

-3 = ^''^ 



(2) 



(3) 



(4) 



and in Eq. ^ each of the P('^(g);/ii) is a Zernike projection 
matrix relating the Zernike decomposition coefficients of the 

wavefront defined on the metapupil (j^*'') and b ^^^^ ^ (gy, hij 
are the Zernike decomposition coefficients of the wavefront de- 
fined on the intersection of the metapupil and the propagation 
path at the Zth layer when looking toward "^(g): 



^lg)(t(g);/!i) = P(t(g);/ii)^ 



(1) 



(5) 



The existence of such P(^(gy,h{) matrices is demonstrated 
inlRa gazzoni et alJ ll 19991) a nd an analytical expression for their 
computation is provided in iFemenia & Brus al ll2003h and de- 
pends exclusively on the geometry of the problem (i.e. on the 
metapupil radius, radius of the circular region defined by the 
intersection of the metapupil and the propagation path and the 
relative position of the circular region with respect to the meta- 
pupil center). The required j"^*^'*!^^ ^ coefficients may be ob- 
tained from Eq. Q by a proper inversion method (e.g. SVD 
or least square estimator, optimal reconstructor or maximum- 
a-priori reconstructor). 



2.2. Implementation in the simulation 

The previous description of modal tomography does not give 
any indication as to how the Zernike polynomial coefficients 



3. Description of the software employed: CAOS 

The main software tool used is the IDL-based Code fo r 
Adaptive Optics Systems (CAOS.'Carbi llet et all2 0Qlb'"l999l). 
The main feature of CAOS is its modular structure: each ele- 
mentary physical process is modeled by a specific set of rou- 
tines (i.e. a module) , e.g. the simulation of atmospheric tur- 
bulence or the propagation of the wavefront through turbulent 
layers. This modular design makes CAOS a very versatile tool 
as it can be extended by the user with new modules as required. 
The other main feature of the CAOS s oftware is its grap hic user 
interface, the "Application Builder" dpini et alJl200fll) . which 
helps to design simulation projects in a very straightforward 
way even for users with no knowledge of IDL. After the graph- 
ical design, the Application Builder automatically generates an 
IDL code which can be run as the desired simulation. In or- 
der to simulate MCAO systems new modules have been devel- 
oped in addition to the modules in the standard CAOS library. 
A dedicated module conducts the calibration of the MCAO sys- 
tem yielding the interaction matrix: at each layer to be recon- 
structed we place a single polynomial Zj at a time and perform 
the geometrical propagation (i.e near-field approximation) of 
the wavefront error corresponding to this Zj along the line of 
sight of each Guide Star (GS) down to the SH sensors. The vec- 
tor of X- and y-slopes from all the SH sensors gives us the col- 
umn of the matrix associated with the j-th Zernike polynomial 
on the considered layer; this process is then repeated for all the 
modes to be reconstructed on each layer The reconstruction 
process is done with a dedicated module which uses a Least 
Square Estimator (LSE) and by inverting with Singular Value 
Decomposition (SVD) the linear relation - IVT^ between 
the measured slopes fi?) and the coefficients (^ ) on each recon- 
structed metapupil. The multi-mirror reconstructor is simulated 
assuming perfect Zernike DMs in the sense that any mirror is 
able to generate a pure Zernike mode up to a given mode j and 
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allows different degrees of correction to be used at different 
layers. 

4. Simulation description & simplifications 

4.1. Generation of Guide Stars (GS) 

NGS are simulated as point-like sources at an infinite distance 
from the telescope and sampling a cylindrical turbulence vol- 
ume. LGS are simulated as extended sources at a finite dis- 
tance so that they sample a cone shaped turbulence volume. The 
width of the LGS is evaluated taking into account the broaden- 
ing due to diffraction by a launching telescope of 0.5 m diam- 
eter, the upward propagation of the laser beam, the perspective 
elongation of the LGS beacon as seen from the observing tele- 
scope and the geometry of the beacon size averaged over the 
range gate (for Rayleigh LGS) or width o f the Na-layer (for Na 
LGS) using the expression presented in lEUerbroek & RigautI 
1I2OOII) : 

where A is the wavelength of observation (A - 589.3 nm for 
Na LGS, for Rayleigh LGS we consider A = 355 nm), dj 
the launching telescope diameter, ro the Fried parameter at the 
WFS operation wavelength, L the distance between the launch- 
ing telescope and the SHS subaperture, h the altitude of the 
LGS and 6h the "range gate". By "range gate" we mean the 
vertical distance of the LGS on which we integrate while it is 
propagating upward. Such a value can be fixed by us when con- 
sidering Rayleigh LGS (by means of e.g. a Pockels cell). In the 
case of a Na LGS it is not possible to fix a range gate since the 
Sodium layer drifts in altitude and thus in the above expression 
we must use the full width of the Sodium layer weighted by 
the fact that the sodium layer has a Gaussian density profile so 
that an appropriate value in the above expression isdh = 5 Km. 
Then the LGSs are projected onto the WFS so that they look 
like circular Gaussian spots of width equal to 6^. This neglects 
the differential elongation due to the finite distance from the 
telescope to the LGS but on the other hand such an approach is 
pessimistic since the LGS is not circular but elongated and in 
practice this involves a loss of resolution only along the elon- 
gation direction. 

4.2. Generation of atmospheric turbulence 

For the purpose of our simulations we have constructed a 7- 
layer profile summarized in Tabled T his profile was built using 
the m easurements of balloon flights dVernin & Munoz-TunonI 
launched from the Observatorio del Roque de los 
Muchachos (ORM). Data from the 6 available balloon flights 
were used to determine an averaged profile in the following 
way; first all profiles from the balloon measurements were nor- 
malized to have the same total ro as the annual mean (i.e. 
ro(500/vm) = 0.15 m) at ORM and th en a continuous averag e 
profile was generated. Then, following Ellerbroek et al. (1994), 
a 7-layer turbulence profile was obtained by performing a dis- 



Table 1. The ORM average balloon 7-layer profile. 
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8 08 X 10" 
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2.15 
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-14 
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1.40 X 10" 


-15 
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4.13 


24.4 



" Height refers to altitude above the telescope. 



Crete fit to the continuous average profile that matches the first 
14 moments of the continuous turbulence profile. 

Atmospheric turbulent phase screen generation was per- 
formed by summing a finite set of Zernike polynomials 
weighted by a set of coefficients drawn from a Gaussian mul- 
tivariate distribution of zero mean and covariance matrix as 
given in .Noll ( 1976), i.e. assuming the turbulence follows 
Kolmogorov statistics. The truncation of the infinite series im- 
plies missing high spatial frequency components. This is easily 
obviated by truncating at a large enough mode y'max- The total 
Strehl ratio due to all missing modes on aU layers is estimated 
to amount to about 10% and this reduction is applied to the 
iso-Strehl maps. A marginal effect is due to the actual statistics 
followed by the Zernike coefficients on portions of the phase 
screens the same size as the telescope pupil. This problem can 
also be circumvented by the use of a large enough number of 
Zernike modes when generating the turbulence phase screens. 

4.3. Propagation through turbulence and wavefront 
sensing 

Propagation of light through atmospheric turbulence is per- 
formed by adding linearly the phase fluctuation produced by 
each turbulent layer which assumes the near-field approxima- 
tion ( Roddier . 1981) and corresponds to neglecting diffraction 
effects in the propagation of the phase of the complex field 
between layers and from layers to the ground (a similar ap- 
proximation is implicit when assuming thin layers in Sect. 14.21 
diffraction effe cts between the to p and the bottom of the layer 
are neglected). Ellerbroek (2001") presents results using a code 
implementing wave-optics for phase propagation and makes 
a comparison of the results obtained under the near-field ap- 
proximation for the case of an MCAO system at the Gemini 
telescope. He concludes that including Fresnel-diffraction in 
the propagation of light involves a phase rms of only a few 
nm (i20 nm) with respect to the results obtained with the ge- 
ometric propagation. Another effect related to not including 
Fresnel diffraction is that we neglect the effect of scintilla- 
tion on wavefront sensing. CAOS implements a wave-optics 
code to simulate the WFS. However the problems of scintil- 
lation are not handled since they do not appear in the near- 
field approximation. Our simulations assume Shack-Hartmann 
WFS for which the performance dependence on the scintilla- 
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tion has been studied in lVoitsekhovich et alJ ( I2OOII) in the case 
of open -loop measurements . A dir ect extrapolation of the re- 
sults of IVoitsekhovich et alJ (1200 ll) indicates that for our case 
with ro(500 nm) - 0.15 m scintillation effects in the WFS will 
involve errors < 8% in the estimation of the high order recon- 
structed modes. A much more important simplification is that 
of no noise in the WFS procedure as we are not concerned at 
this stage with problems associated with sky coverage and by 
not considering noise we are in the limit of very bright objects. 

4.4. Least Square Estimator (LSE) reconstructor 



SVll) Threshold Sckiction 



Several works (IFusco et"al] Il99^ Brusa et all 
iFhcker et alJ l2000t IFusco et alJ 1200 it 



pooot 

Conan et alJ l200lh 

have focused on estimators for MCAO which assume a real- 
time knowledge of the statistics of phase perturbations and the 
turbulence vertical distribution to derive optimal reconstructors 
operating in open-loop: the so-called Maximum a posteriori 
(MAP) estimators. In conditi ons of low Signal-to-Noise Ratio 
(SNR) it has been shown jFlicker et all l2000t IConan et al] 
[2001 ; LeLouarn 20o2) that these optimal reconstructors 
perform better than the SVD-based reconstruction while 
yielding similar results in the case of high SNR values. Since 
we work in the limit of very bright objects (see Sect. 14. 3> 
nothing justifies using MAP estimators here except for the 
advantage of including knowledge of the vertical distribution 
of turbulence. On the other hand, it is clear that a MCAO 
system will be operated in closed-loop which imposes no 
problem for the implementation of the LSE reconstructor 
while, to our knowledge, statistics of the residual signal on the 
WFS have not been derived analytically to be implemented in 
the design of the MAP reconstructor. 



4.5. Perfect Zernike defer mable mirrors 

The DMs considered in these simulations are assumed to be 
able to reproduce exactly the Zernike modes being recon- 
structed, that is, if the command is that of generating ajZj (with 
j = 2, . . . , 153) on DMi then this is exactly what we get (idem 
with DM2 but j -2, . . . ,91) . This assumption neglects the in- 
fluence functions of real DMs which in addition to the com- 
manded Zj mode will probably introduce other modes at a much 
lower level. Realistic influence functions will be considered in 
future work using more detailed numerical simulations. 

4.6. Open-loop operation. 

Our simulation assumes open-loop operation. In doing this we 
can easily perform a statistical analysis by generating inde- 
pendent input turbulent wavefronts and retrieving the response 
of the MCAO system. Open-loop operation neglects degrada- 
tion of the MCAO system performance due to time-delay er- 
rors while on the other hand it gives an underestimation of the 
performance of the system as compared to closed-loop oper- 
ation since the WFS is more affected by problems due to the 
restricted FoV in the SHS and the extended size of the GS. 



0.50 



a 0.45 



0,40 



0.35 



0.15 



0.10 



a 0.05 



0.00 





— l.OE-4 




■-0 1.5E-4 


— 


2.8E-4 


— 


- O5.0E-4 


0----- 


-0 l.OE-3 




2.3E-3 







— 1 


OE 


-4 





1 


dE 


-4 




--^2 


8E 


-4 


0--- 


---0 5 


OE 


-4 




-0 1 


OE 


-3 


<y - 


-^2 


3E 


-3 




15 20 25 30 
Radius (arcsec) 



Fig. 1. Example of SVD threshold selection for the case of 
3NGS on a circumference of 45" with DM2 conjugated to 
4.7 Km.The normalized S rms curve corresponding to thr2 - 
1.5 X 10"** is slightly displaced upward so as not to be super- 
imposed on the curve corresponding to thrj - 2.8 x 10"^. Here 
we selected the thrj, curve as the one with the best compromise 
between the highest S ratio value and the lowest S rms. 



5. Simulation results 

In the case of single-conjugate AO with moderate-good cor- 
rection, the performance may be characterised adequately by 
a single parameter, namely the on-axis Strehl ratio. The Strehl 
ratio can be maximised by choosing a threshold value in the 
SVD reconstruction; singular values lower that the threshold 
are set to zero and this eliminates modes to which the system is 
insensitive. In the case of an MCAO system, the performance 
to be determined is the Strehl ratio over the field of view. We 
find in our simulations that the Strehl ratio maps are sensitive 
to the SVD threshold (for example, see Fig.|5J- Since the SVD 
threshold cannot be selected a priori we need criteria for select- 
ing the optimal threshold values for each configuration. In or- 
der to characterize SR maps, we extract two parameters which 
quantify the average Strehl ratio over the field and the amount 
of variation of the Strehl ratio over the field. For a given Strehl 
map we determine the value of the Strehl ratio averaged on con- 
centric annular radii (the first annulus is the inner circle of ra- 
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dius 7" and concentric annuli have an outer radius 2.5" larger 
than the inner radius) and the root mean square deviation of 
Strehl in each annulus with respect to the averaged Strehl ra- 
tio on the circle of radius equal to the annulus outer radius. 
As shown in the example in Fig. ^ we obtain two families of 
curves: the first set of curves shows for each SVD threshold 
value the evolution of Strehl averaged on annuli as a function 
of the annulus outer radius while the second family of curves 
shows the evolution of Strehl nns. This process is conducted 
on contour plots such as those depicted in Fig.|2l To produce 
these contour plots we interpolate the reconstructed FoV sam- 
pled at 15x15 locations in the simulation to a finer grid. Due to 
this interpolation when we get close to the edges we observe a 
large gradient in S. To avoid a significant bias in the SVD se- 
lection process we generate curves of averaged Strehl ratio and 
Strehl nns on concentric annuli up to a maximum outer radius 
of about 85% of the whole reconstructed FoV radius. Once we 
have generated the above sets of curves for each GS configu- 
ration and DM2 conjugation altitude the SVD threshold value 
selection is based on (i) having flat curves of averaged Strehl 
within concentric annuli (in general we observe curves which 
are mostly monotonically increasing or decreasing, so that we 
select those curves with the lowest peak-to-valley values and/or 
the lowest gradient on annulus averaged Strehl), (ii) having the 
smallest Strehl nns and (iii) in those situations where several 
threshold values yield similar results based on the previous two 
points we select the largest SVD threshold value. This is justi- 
fied in terms of SNR considerations. 

5.1. Performance using Natural Guide Stars 

Fig. |5]shows K-band Strehl maps for the case of having three 
natural guide stars evenly spaced on a circumference of radius 
30" in order to reconstruct a circular field of view of 1' diam- 
eter for several SVD threshold values. The panels correspond 
to setting different singular value thresholds; it can be seen that 
the performance is fairly constant until a large number (15) of 
SVD modes are rejected, some of them with important con- 
tributions of absolute Zernike tilts on each of the two recon- 
structed layers. When only a few SVD modes are rejected (i.e. 
from THR=7 x 10"^ to THR=3 x 10 involving 2 to 7 SVD 
modes rejected) the Strehl ratio peaks at the positions of the 
guide stars as would be expected considering the least-squares 
reconstruction process employed while the Strehl ratio across 
the reconstructed FoV is relatively uniform when rejecting 8 to 
10 SVD modes (THR=4 x 10 ^* to THR^IO^^). This is a gen- 
eral behavior observed in our simulations: uniformity of the 
Strehl ratio across the reconstructed FoV is achieved at the cost 
of lowering the peak Strehl ratio at the GS positions by an ad- 
equate SVD filtering. Fig.|3]shows the best SVD Strehl maps 
when the three natural guide stars are on a circumference of 
radius of 45", corresponding to a reconstructed FoV of 1.5' di- 
ameter, and DM2 conjugated to 6.1 and 10 Km above the tele- 
scope pupil. We observe that with DM2 conjugated to 10 Km 
the average Strehl ratio is significantly smaller than when DM2 
is conjugated to 6.1 Km while the variation of the Strehl ra- 
tio across the field is more pronounced for DM2 conjugated to 



10 Km. Making measurements on the maps it is found that the 
average Strehl ratio drops from 0.43 to 0.20 and the normalised 
Strehl nns increases from 16% to 48% (see also Fig.l^. The 
Strehl nns value when DM2 is conjugated to 10 Km (and to 
13 Km) is due to such a strong center-to-edge Strehl variation 
that there is little sense in talking of an average reconstructed 
Strehl ratio for this scenario. 

By varying the DM2 conjugation altitude in our simulations 
we obtained that for constellations of 3 NGS the optimal conju- 
gation altitudes are approximately 8.5 Km and 6 Km for recon- 
structed FoVs of r diameter and 1.5' diameter, respectively. 
These optimal altitudes are surprisingly low when we consider 
the vertical distribution of turbulence in the 7-layer model in 
Table [l] We would expect the pupil-conjugate DM to correct 
the turbulence in the four lower layers (up to 5 Km), and the 
second DM would therefore be positioned to correct t he turbu- 
lence in the three upper layers (10-18 Km). Tokovinin & Viard 
(2001) published expressions allowing optimal conjugation al- 
titudes to be calculated for up to three mirrors for a given turbu- 
lence profile. Applying these expressions to our case (i.e. two 
DMs, balloon 7-layer model), we obtain an optimal DM2 alti- 
tude of 13 Km. We believe the reason for such a disagreement 
has to do with the strong assumptions in the aforesaid theo- 
retical work: the authors assume that whatever the altitudes of 
conjugation for the DMs, tomography is able to perfectly re- 
construct the wavefront at those altitudes so that the only effect 
considered is the correction of a continuous turbulence profile 
by a finite number of DMs. In fact, the further from the tele- 
scope we conjugate DM2, the less overlap between the GS foot- 
prints at the metapupil at the DM2 conjugation altitude. This ef- 
fect will tend to push down the optimal conjugation altitude for 
DM2. This statement is supported by the observation that the 
optimal conjugation altitude is lower for the case of a FoV of 
1.5' than for the case of a FoV of 1'. In order to further test this 
hypothesis, we carried out simulations using five natural guide 
stars and keeping the rest of the parameters of the simulation 
as in the case of 3 NGS. Increasing the number of NGS in our 
constellation from 3 to 5 has the effect of increasing the overlap 
between their footprints and thus we decrease the importance of 
voids in the metapupil for the wavefront reconstruction at that 
range. A comparison between the results achieved with constel- 
lations of 3 NGS and constellations of 5 NGS is given in the left 
top and left bottom panels in Fig. |2 where we plot the average 
Strehl ratio on the reconstructed FoV as a function of the DM2 
conjugation altitude and the Strehl ratio nns variation, respec- 
tively (the right top and right bottom panels of Fig. |4] show the 
results obtained when considering different LGS constellations 
plus an on-axis NGS and which will be the subject of the fol- 
lowing sections). We can see in Fig.0]that when increasing the 
number of NGSs from 3 to 5 NGS the conjugation altitude of 
DM2 also increases and this effect can only be due to a better 
sampling (and reconstruction) of the turbulence at the height at 
which DM2 is conjugate d. From these resul t s we c onclude that 
the effect neglected by Tokovinin & Viar J 1200?) is very im- 
portant when determining the best DM2 conjugation altitude. 

We also compared two different geometries employing 5 
NGS: placing all NGS on the vertices of a pentagon or putting 
4 of them on the corners of a square and one at the center. It 
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Fig. 3. Strehl maps for a reconstructed 1.5' FoV for DM2 conjugated to 6.1 Km and 10. Km for a constellation of 3NGS on a 
circumference of 45" radius. These maps correspond to the best SVD case for each DM2 conjugation altitude. 
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Fig. 4. Summary of simulation results. Top panels show the dependence of the MCAO average Strehl ratio on the conjugation 
altitude of DM2 ■ The bottom panels give the normalized nns of the averaged Strehl ratio on the 83% radius annulus (see main 
text in Sect.|5lfor further details). 



turns out that it is the latter configuration which yields better 
results (see top left panel in Fig.|?]and iso-Strehl contour plots 
in Fig. |5} despite the fact that it leaves unsampled a slightly 
bigger portion of the metapupil. This suggests that not only 
the fractional sampling of the metapupil is important but also 
the distance of points within the reconstructed FoV to the GS. 
Although we have not conducted any test to check our expla- 
nation of this effect we believe that this may be due to the LSE 
reconstruction in which by definition the best reconstructed po- 
sitions are those corresponding to the GS positions while for 
a MAP reconstructor we may encounter a different behavior 
since it is built by optimizing the reconstruction at a set of po- 
sitions which do not have to be coincident with those of the 
GS. 

5.2. Performance using Sodium Laser Guide Stars 

The probability of finding constellations of 3 (not to mention 
5) bright natural guide stars arranged to nicely span fields of 
1 - 1 .5 ' is very low. As an indication, using the Guide Star 
Catalogue-ll' we estimate that the probability of finding three 

' Copyright Association of Universities for Research in 
Astronomy, Inc. Further details and access to database at 
|http://www-gsss. stsci.edu/gsc/gsc2/gsc22_release_notes.htm_ 



stars brighter than wr = 16, which is an optimistic value for 
limiting guide star magnitude, is approximately 3% for a field 
of 2' at galactic coordinates (/, b) = (90° , 40°), and this value 
will be further reduced when requiring the 3 NGS to be prop- 
erly arranged. In view of these considerations it will be nec- 
essary to employ laser guide stars (LGS) in order to provide 
reasonable sky coverage for star-oriented MCAO systems, at 
least on 8-lOm telescopes. 

It has already been pointe d out by different authors (see 
e.g. Eller broek & Rigauti200ll) that the indetermination of the 
tip-tilt of individual LGS leads to tip-tilt anisoplanatism in an 
LGS-based MCAO system. The tip-tilt anisoplanatism may be 
solved in a natural way by detec ting the tip-tilt of several NGS 
in the FoV ('Rigaut et al. 2000V In the Gemini-South MCAO 
system it is planned to reconstruct the tip-tilt using three dim 
(i.e. mR ~ 19 - 20) natura l guide stars and th e sky coverage 
in the K band is predicted jRigaut et al ''200Cf) to be 15% and 
80% at Galactic pole and b - 30°, respectively. Alternative 
schemes can be derived by considering that tip-tilt anisopla- 
natism mainly arises from turbulence-induced quadratic wave- 
front eiTors (i.e defocus and astigmatism) at altitudes which 
cannot be determined from the multiple-LGS measurements. If 
quadratic measurements are made simultaneously with guide 
stars at different ranges, then it may be possible to correct the 
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Fig. 5. Checking two constellation geometries employing 5 NGS. On the left panel the Strehl map when the 5NGS are evenly 
arranged on the vertices of a regular pentagon within a 1.5' FoV. On the right panel the results when 4 NGSs are placed on the 
vertices of a square and 1 NGS is placed on the FoV center. Both panels correspond to having DM2 conjugated to 10 Km above 
the telescope pupil. 



1' FoV using 4 LGS 1.5' FoV using 5 LGS 




-40 -20 20 40 -40 -20 20 40 



Fig. 6. Left panel shows the best SVD Strehl map in a reconstructed FoV of 1' diameter with DM2 conjugated to 10 Km and 4 
LGS on a circumference of 41.46" radius plus one on-axis bright NGS (i.e. wave-front sensing on the on-axis NGS is up to the 
same order as with any LGS). Right panel is the Strehl map corresponding to the reconstruction of a 1.5' diameter FoV with 
DM2 conjugated to 10 Km and 5 LGS on a circumference of 56.46" radius plus one on-axis bright NGS. 



tip-tilt anisoplanatism. One approach along these lines is to 
make wavefront measurements on a natural guide star simul- 
taneously with the laser guide stars. We determine the perfor- 
mance of such an approach here for fields of view of diameter 
r and 1.5'. In the case of the smaller field of view we employ 



four Na LGS on a circumference spanning the field of view 
plus one on-axis NGS. In the case of the 1.5' field of view we 
consider two Na LGS configurations; (i) four LGS spanning 
the field of view plus one on-axis NGS (ii) five LGS on a cir- 
cumference spanning the field of view plus one on-axis NGS. 
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The circumferences on which the LGS are placed are some- 
what larger than the fields of view so that the metapupil sizes 
are the same when using cone shaped beams from the LGS as 
when using cylindrical-shaped beams from NGS placed at the 
edge of the field of view. Fig.|6]shows the best Strehl maps for 
the 1' FoV (left panel) and 1.5' FoV (right panel), employing 
constellations of 4 ofF-axis LGS+ on-axis NGS and 5 ofF-axis 
LGS -I- on-axis NGS, respectively . For the 1' FoV case the av- 
erage Strehl ratio is similar to that obtained on the same field of 
view using 3 NGS although it exhibits a larger peak-to-valley 
value. The rms Strehl variation is similar to the 3 NGS case but 
it does not rise steeply when DM2 is conjugate to an altitude 
higher than the optimal which for the LGS case is approxi- 
mately 10 Km while in the case of NGS constellation it was 
about 8.5 Km. In the case of the 1 .5' FoV the average Strehl ra- 
tio is almost the same for both 4LGS and 5LGS configurations. 
The optimal altitude occurs in the range 8-10 Km. If DM2 is 
conjugate to an altitude higher than the optimal then the aver- 
age Strehl ratio is higher and the rms variation is lower when 5 
LGS are employed rather than 4 LGS. 



5.3. Hybrid Rayleigh-Na system performance 

The natural guide star employed to measure quadratic modes 
has to be bright, thereby again compromising sky coverage. 
An alternative approach to solving the tip-tilt anisoplanatism 
caused by the quadratic modes consists of employing both 
Rayleigh and Na LGS as well as a single nat ural guide star 
which is onl y used to determine global tip-tilt (Fe menfa et alJ 
I2OOL EUerbroek & Rigaut 200 L Femema 2003) . Here we ap- 
ply the algorithm described in Femenial ]2003l) whose starting 
point is to realize that unresolved second order modes in LGS- 
based tomography can be reconstructed using an LGS con- 
stellation made at several altitudes (e.g. Rayleigh + Sodium 
LGSs, or Rayleigh LGSs at different altitudes) as shown by 
iBrusa et aL C2000J . The reconstruction process consists of two 



stages in which barycenter tip-tilt is measured from a NGS and 
in parallel we perform tomography with the LGS constella- 
tion placed at at least two different altitudes. From the LGS 
tomography we obtain a reconstruction of all modes starting 
from the Zernike second-order modes and this information is 
used together with the NGS barycenter tip-tilt to obtain an esti- 
mate of the global Zernike tilts in the entire reconstructed FoV. 
The hybrid technique may be implemented in different ways, in 
particular Femenfa (200,3) considers a single on-axis Rayleigh 

LGS on which to conduct h igh-order wave-front sensing while 
1 1 1 

EUerbroek & Rigaut (2001) consider several off-axis Rayleigh 
LGS coupled to low-order WFS. Other parameters to be opti- 
mized include the height of the Rayleigh LGS and the wave- 
length for the Rayleigh LGS. Given the number of new pa- 
rameters to be investigated the results we present here are by 
no means exhaustive but only preliminary to give an idea of 
the potential of the technique. We assume two configurations 
of 4 Na LGS and 5 Na LGS evenly arranged on a circumfer- 
ence of 56.45" radius (so that the reconstructed FoV is 1.5') 
plus an on-axis Rayleigh LGS. The Na LGSs are observed at 
589 nm and the Rayleigh LGS at 355 nm with 16 x 16 sub- 



aperture SHSs. The rest of the simulation parameters (atmo- 
spheric profile, number of CCD pixels and SHS FoV, no noise 
in the SHS, number of reconstructed modes, etc.) are the same 
as in the simulations with NGS and LGS+on-axis NGS. The 
top panels a) & b) in Fig.0shows the reconstruction across the 
1 .5' FoV for the cases of 4 Na LGS and 5 Na LGS, respectively. 

On comparing the results from the constellation using 4 Na 
LGS against the results using 5 Na LGS we see a limited gain. 
This limited performance gain is a consequence of the limited 
reconstruction along the on-axis direction which in turn is used 
to reconstruct the Zernike tip-tilt modes across the entire re- 
constructed FoV. This statement is supported by the levels of 
reconstruction attained if we assume perfect reconstruction of 
the Zernike tip-tilt modes across the entire FoV as depicted in 
panel c) in Fig.0for the case of 4 Na LGS at 56" and in panel 
d) in Fig. 0for the case of 5 Na LGS at 56"; in both cases 
also considering on-axis Rayleigh LGS at 30 Km and on-axis 
Tip-tilt NGS. When comparing in Fig.0the bottom panels (per- 
fect correction of Zernike tip-tilt across the entire FoV) against 
the top panels (actual correction of Zernike tilts) we verify that 
in the latter figures there is a small amount of tilt anisopla- 
natism caused by a worse coiTection of higher-order modes (es- 
pecially second-order modes) in the on-axis direction. This in 
turn translates into tilt anisoplanatism when using those higher- 
order modes to reconstruct off-axis Zernike tilts. This situation 
also appeared when dealing with 5 NGS where we observed 
that the best configuration is that corresponding to 4 NGS on 
the corners of a square plus one on-axis NGS rather than all 5 
NGS on the vertices of a pentagon. Based on this conjecture it 
would be more effective having 4 Na LGS distributed on the 
vertices of a square plus one on-axis Na LGS than having all 5 
Na LGS distributed on the vertices of a pentagon. 

6. Choosing the DM2 conjugation altitude 

As seen in Fig.|4]the optimal DM2Conjugation altitude depends 
very strongly on the GS nature and configuration. In the case of 
observing 3NGS on the vertices of equilateral triangles on a cir- 
cumference of radius 30" or 45" the optimal altitudes are about 
8.5 and 6.0 Km, respectively. When considering 5 NGSs the 
optimal DM2 conjugation altitude is lOKm. When hybrid LGS 
constellations and a single tip-tilt NGS are used aiming at a re- 
constructed FoV of 1.5' the best DM2 conjugation altitude oc- 
curs in the range [8,1 l]Km with little variation of performance 
within that range leaving the final decision to opto-mechanical 
considerations. In order to find an optimal DM2 conjugation al- 
titude that best suits all possible GS configurations cases in this 
work we consider the following figure of merit: 



/(ft) 



Z 



S max.i S \{K} 



' max.i 



where h denotes a generic DM2 conjugation altitude, the index 
i = 1 , . . . , 9 runs over all the GS configurations in this work, 
Wi weights different GS configurations, ^max.i is the maximum 
Strehl ratio achieved with the ith GS configuration and S i{h) is 
the achieved Strehl ratio with the /th GS configuration when 
conjugating DM2 to altitude h above the telescope entrance 
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Fig. 7. Strehl maps for two LGS geometries employing the hybrid LGS constellations. Panel a) coiTesponds to using 4 Na LGS 
at 90 Km with an ofF-axis angle of 56.46" (in order to provide a coiTected 1.5' FoV) and an on-axis Tip-Tilt NGS. Panel b) 
corresponds to using 5 Na LGS at 90 Km with an ofF-axis angle of 56.46" plus an on-axis Rayleigh LGS at 30 Km and an on-axis 
Tip-Tilt NGS. Both cases consider DM2 conjugated to 10 Km. Panels c) and d) are the the same as panels a) and b), respectively 
except that tip-tilt is perfectly removed. This ideal exercise supports our interpretation on the limited level of reconstruction being 
mainly caused by residual tip-tilt modes. 



pupil. As the mimization process will depend on {wi} we 
have considered the sets of weights listed in Table |2l the first 
weighting set (i.e. column wiin Table |2j gives equal impor- 
tance to all GS configurations, weighting set {W2} is focused on 
avoiding system complexity and thus assigns higher weight to 
simpler systems while weighting sets {wy,} and {W4} are derived 
from simple sky coverage considerations based on star counts 
from the Guide Star Catalogue-2 at (/, b) = (90° ,40°) and as- 
suming that for the NGS we require at least < 16 when 
it is used for high-order mode reconstruction and wr < 18.5 
when it is only used for tip-tilt reconstruction (i.e. as in the 
case of the hybrid laser configurations). The difference between 
weighting sets {W3} and {W4} is that the former is proportional 
to the probability of finding the required number of NGSs in the 
GS configuration and the second is proportional to the square 
root of that probability (i.e. W44 = yw3,i where index / runs 
over the different GS configurations under study). This choice 



for {w3)and {W4} is very qualitative but given the small differ- 
ence between the results obtained with them (see Table |3} we 
are confident that a detailed sky coverage analysis is likely to 
differ little from the results obtained here. Thex^W minimiza- 
tion process yields an "optimal" DM2 conjugation altitude as a 
function of the weighting set. The first and second columns in 
Table|3lgive the DM2 conjugation altitude for which a given GS 
configuration attains its maximum Strehl ratio and the Strehl 
ratio itself. Columns 3 to 6 give Strehl ratios for each GS con- 
figuration when the DM2 is conjugated to the joint optimal con- 
jugation altitude according to the weighing set under consider- 
ation (i.e 5opt,iCorresponds to using wi); the value of such opti- 
mal conjugation altitude for each weighting set being displayed 
in parenthesis below each column title. The last two columns 
in Table |3give the Strehl ratio and fractional decrease with re- 
spect to Smax, respectively, when DM2 is conjugated to 8 Km. 
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Table 3. Finding the optimal DM2 conjugation altitude. 
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Table 2. Different normalized-to-maximum weighting sets 
used to obtain the DMaconjugation altitude. Values marked 
with superindex " should be further reduced when requmng 
proper arrangement of the NGSs. 
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7. Conclusions 

Modal MCAO in the "classical" star-oriented fashion holds the 
promise to reconstruct FoVs as large as 1 .5' in diameter. Ideally 
one would like to rely on NGS constellations but the chances 
of finding a suitable one is extremely low even for the sim- 
plest case of 3NGS. Thus if one wants to build a star-oriented 
MCAO system which can be used over a significant fraction 
of the sky one must necessarily consider LGS-based MCAO 
systems. However, one then has to face the tip-tilt problem 
and the indetermination of second-order modes. The simplest 
approach to LGS-based MCAO systems is that of relying on 
a bright enough NGS on which one could perform not only 
tip-tilt & second-order sensing but also high-order WPS. This 
yields poor results and the initially surprising result that per- 
formance does not improve significantly when increasing the 
number of Na LGSs from 4 to 5. In any case, whether we use 4 
Na LGS or 5 Na LGS the reconstructed FoV is strongly affected 
by tip-tilt anisoplanatism and the reconstruction is extremely 
sensitive to the SVD filtering process which in practice wiU be 
determined by the SNR in the WFS measurements. We also re- 
mark that such a naive approach also has small sky coverage 
as it requires a star which is bright enough to obtain high-order 
wavefront measurements. However, with hybrid LGS constel- 
lations (Sect. 15.31 it may be possible to obtain good levels of re- 
construction. Taking the sky coverage values obtained with the 
three tip-tilt NGS scheme considered by Gemini (see Sect. 15. 2t 
and assuming equivalent magnitude requirements for the tip- 



tilt NGS with hybrid LGS constellations it is expected that this 
approach can be applied to a larger fraction of the sky as only 
a single tip-tilt NGS is required. There is a lot of room for pos- 
sible improvements to the technique used ( lFememal2003l) . 

The optimal DM2 conjugation altitude depends strongly 
on the guide star constellation used. When considering NGS 
constellations our MCAO simulati ons give results in di s agree- 
ment with the theoretical work by 'Tokovin in & Viar3 ( I2OOII) 
which would place D Mj at about 1 3 Km above telescope pupil 
( iDevanev et alJ2002h . The reason for such disagreement is be- 
lieved to be caused by the strong assumptions in these theo- 
retical works resulting in overestimating the DM2 conjugation 
altitude. From our study based on detailed numerical simula- 
tions we consider that a suitable DM2 conjugation altitude for 
a 10-m class telescope is around 8 Km and thus we are con- 
sidering such an option for the design of the AO system at the 
GTC telescope which considers an initial AO system with a 
single corrector but upgradeable to a dual-conjugate system. 
It has been argued that the aforementioned theoretical works 
would be applicable to the case of extremely large telescopes 
(30-100 m diameter) as the main limiting assumption (i.e. full 
metapupil coverage by overlapping the GS beams at the recon- 
structed layers) is relaxed and it should be possible to fulfill the 
requirement of perfect wavefront reconstruction at each meta- 
pupil. However, one should consider that the larger the tele- 
scope diameter, the larger the number of modes required in 
order to reach the same Strehl Ratio as with a smaller tele- 
scope diameter. It follows from this consideration that with ex- 
tremely large telescopes, although the metapupil coverage is 
nearly complete, the importance of small gaps in the metapupil 
coverage is much more important, this implying that the as- 
sumptions in the above theoretical works may not be realistic 
even for extremely large telescopes. 

Our last conclusion regards a rule of thumb indicated by 
our simulations: given a fixed number of GS (either LGS or 
NGS) the best MCAO performance occurs for the constellation 
configuration that minimizes the distance from any point in the 
FoV to a GS position. As we have seen although 5 GS on the 
vertices of a pentagon give a larger covered surface of the meta- 
pupil than 4 GS on the square vertices plus one on-axis GS, the 
latter configuration yields the better reconstruction. 

We would like to end with a comment on the turbulence 
profile used in our simulations. It has been derived from only 
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6 turbulence profiles at the ORM obtained with balloon flights 
which is a rather small sample of turbulence profiles. A larger 
set of turbulence profiles would be desirable in order to obtain 
a more significant average profile at ORM. In any case, all the 
simulations were also conducted on a modified version of the 
average profile and the optimal conjugation altitude of DM2 
remained essentially the same. 
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